Astronomy & Astrophysics manuscript no. aruiz SEDhlirgs 


©ESO 2010 


March 3, 2010 





Spectral Energy Distribution of Hyper-Luminous Infrared Galaxies 

A. Ruiz^'^, G. Miniutti^, F. Panessa^, and F.J. Carrera' 



Institute de Fi'sica de Cantabria (IFCA), CSIC-UC, Avda. de los Castros, 39005 Santander, Spain 

Istituto Nazionale di Astrofisica (INAF), Osservatorio Astronomico di Brera, via Brera 21, 20121 Milano, Italy 

LAEX, Centra de Astrabiologia (CSIC-INTA); LAEFF, P.O: Box 78, E-28691 Villanueva de la Caiiada, Madrid, Spain 

Istituto Nazionale di Astrofisica (INAF), lASF-Roma, Via Fosso del Cavaliere 100, 1-00133 Rome, Italy 

e-mail: angel.ruiz@brera.inaf.it 



o 

(N 



O 
U 

O 



Received March 31, 2009; accepted March 2, 2010 

ABSTRACT 

Aims. The relationship between star formation and super-massive black hole growth is central to our understanding of galaxy for- 
mation and evolution. Hyper-Luminous Infrared Galaxies (HLIRGs) are unique laboratories to investigate the connection between 
starburst (SB) and Active Galactic Nuclei (AGN), since they exhibit extreme star formation rates, and most of them show evidence of 
harbouring powerful AGN. 

Methods. Our previous X-ray study of a sample of HLIRGs shows that the X-ray emission of most of these sources is dominated by 
AGN activity. To improve our estimate of the relative contribution of the AGN and SB emission to its total bolometric output, we have 
built multi-wavelength (from radio to X-rays) spectral energy distributions (SEDs) for these HLIRGs, and we have fitted standard 
empirical AGN and SB templates to these SEDs. 

Results. In broad terms, most sources are well fitted using this method, and we found AGN and SB contributions similar to those 
obtained by previous studies of HLIRGs. We have classified the HLIRGs SEDs in two groups, named class A and class B. Class A 
HLIRGs show a flat SED from the optical to the infrared energy range. Three out of seven class A sources can be modelled with a pure 
luminosity-dependent QSO template, while the rest of them require a type 1 AGN template and a SB template. The SB component is 
dominant in three out of four class A objects. Class B HLIRGs show SEDs with a prominent and broad IR bump. These sources can 
not trivially be modelled with a combination of pure AGN and pure SB, they require templates of composite objects, suggesting that 
> 50% of their emission comes from stellar formation processes. 

Conclusions. We propose that our sample is actually composed by three different populations: very luminous QSO (class A objects 
with negligible SB contribution), young galaxies going through their maximal star formation period (class A objects with significant 
SB emission) and the high luminosity tail of ULIRG population distribution (class B sources). 
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1 . Introduction quire that most AGNs in the Universe are obscured jUeda et al] 

25 [2003; Gilli et al. 2007), i.e. most of the accretion power in the 

^ During the last decade, the hypothesis that Active Galactic Universe is absorbed an d then re-emitted in the infrared (IR) 

Nuclei (AGN) and galaxy formation and evolution are closely bands (Fabian & Iwasawalll999l) . 

O related has been supported by a growing body of observational ir and X-ray observations are therefore essential to under- 

O evidences. On one hand, most galaxies have been shown to har- stand the phenomena of star formation and AGN, as well as their 

bour a central super-massive black hole (Kormendy & Gebhardt interplay and connection. Fortunately, nowadays we have pow- 

^ ; |2001|) w hose mass is correlated with that of the host g alaxy erful tools to observe the Universe in both energy ranges, like 

. spheroid (|Magorrian et alj|1998t [McLure & Dunlop||2002|) and, Chandra, XMM-Newton, Spitzer, AKARI or Suzaku. Different 

on the other hand, the evolution of cosmic st ar formation and of sti-ategies can be employed to investigate the IR/X-ray syn- 

luminous AGN activity appear rather similar (Franceschinietal. grgy and its effect on the AGN-galaxy co-evolution, e.g. by 
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11999^ .Silverman et al, 2005). These hints clearly suggest a con- multi-wavelength surveys like GOODS, AEGIS or COSMOS 

nection between the growth of the central black hole through jDickinson etlW^lM^ iDavis et al.ll2007t IScoviUe et alJIIO^Th . 

accretion and the growth of the spheroid through star formation, by targeted MIR observations of peculiar X-ray sources like X- 

The observational study of these two phenomena needs pen- ray absorbed broad line QSO ( Stevens et al. 2005; Page et^ 

etrating radiation like X-rays, mid-infrared (MIR), fai--infrared 2007), and by tai-geted X-ray observations of MIR/FIR- 

(FIR) or sub-mm. On one hand, star formation takes place in emit ting objects lik e Ultraluminous Infrared Galaxies (ULIRGS, 

heavily obscured environments. Primary radiation is then repro- Franceschini et alj |2003; Teng et al. 2005]) or Hyperluminous 

cessed by dust and re-emitted in the MIR-FIR band. X-ray emis- infrared Galaxies (HLIRGS, |Wilman et al] |1998|; |Ruiz et al.| 

sion from starburst (SB) activity is enhanced by energetic phe- .20071). 

nomena related to the final stages of stellar evolution, e.g. su - ULIRGs are a family of galaxies with IR luminosity Lir > 

pernova remnants or X-ray binaries ("Persic & Rep haelill2002h . 10'^ Lq, whose bolometric output is dominated by the emis- 

On the other hand. X-ray emission is the signature of AGN sion in the IR waveband (see Lonsdale et al. 2006 for a com- 

activity, produced by black hole (BH) growth through accre- plete review). X-ray and IR data clearly suggest that these 

tion. However, synthesis models of the X-ray background re- sources are powered by SB and, in some cases (~ 50%), by 
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AGN s (iFarrah et alll2003l: iFranceschini et alJl2003t iTeng et alJ 
l2005t iNardini et al.ll2008h The fraction of ULIRGs hosting an 
AGN increa ses with increasing IR luminosity (Veilleux et al. 
I1995L Il999h . Most of these objects are in interacting systems, 
i.e. ULIRGs are most likely triggered by mergers of galaxies 
(IFarrah et al.ll200ll I Veilleux et al.ll200l . 

HLIRGs present an IR luminosity Ljr > 10'^ Lq. These 
are among the most luminous objects in the Universe. Assuming 
that the FIR emission above 50 jum is dominated by SB, their 
estimated star formation rates (SFR) are > 1000 Mq yr"' 
(Rowan-Robinson 2000). IR and o ptical observa t ions support 
that in ost harbour an AGN dVerma et alJ I2002L IFarrah etaP 
I2002ah . although the main power source is still controversial. 
As HLIRGs could represent the most vigorous stage of galaxy 
formation, they are unique laboratories to investigate extremely 
high stellar formation, and its connection to super-massive black 
hole growth. 

Only about a third of HLIRGs are located in interacting sys- 
tems (Farrah et al. 2002b), so a considerable number of these ob- 
jects can not be classified just as objects in the brightest end of 
the ULIRG population. They could be very young ga laxies expe- 
rienc ing their major episode of star formation ( Row an-RobinsonI 
12000 ). or may be a completely new class of objects, e.g. a tran- 
sient IR-lumi nous p hase in quasar evolution (.Farrah et al..2002aL 
IStevens et a02005l) . 

X-rays are a very convenient tool to disentangle the relative 
contribution of SB and AGN to the total output of HLIRGs. Only 
a few of these objects had been st udied i n X-rays (jW ilman et al. 
119981 l2003t llwasawa et aDl2005l) before iRuiz etall (.2007.) pre- 
sented the first systematic study of these sources in the X-ray 
band. 

A sample of 14 HLIR Gs was obse rved by XMM-Newton 
and 10 were detected (Ru iz et al.ll2007l) . All of them show an 
AGN-dominated X-ray spectrum. We find X-ray thermal emis- 
sion associated with SB for just one source, while all ULIRGs 
show a SB component in their X-ray spectra (IFranceschini et all 
1200 3). The much brighter AGN emission probably hides the X- 
rays originated in the SB (if this component actually exists). The 
IR luminosity of most HLIRGs of the sample is consistent with 
an AGN or igin, but it is systematically over that ex pected for a 
local QSO (lElvis et alJll9"94l:lRisahti & Elvisll2004l) of the same 
X-ray luminosity. This IR excess could be due to X-ray obscu- 
ration, SB emission or may be due to an intrinsic difference be- 
tween the spectral energy distribution (SED) of AGN in HLIRGs 
and the SED of local QSO. 

To clarify these questions a proper study of the SED 
of these objects is nee ded. Several studi e s of H LIRGs SED 
have been published ( Rowan-RobinsonI l2000t [Verma et_al.. 
l2002t IFarrah et al.) l2002al) . but thev were always limited to 
the IR energy range. These studies apply a two compo- 
nent model (AGN-hSB) to reproduce the IR emission, us- 
ing r adiative transfer models (RTM) fo r the AGN dust 
torus dEfstathiou & Rowan-Rob insonI fl995^ i Rowan-RobinsonI 
1995!) and the SB (Efstathiou et al. 200^) components. 
Rowan-RobinsonI (12000 ) studied a sample of 45 HLIRGs, find- 
ing a continuum distribution in the relative contribution of 
the AGN and SB components, from pure starburst to pure 
AGN, with most ob jects being composite. On the other hand, 
IFarrah etalJ(l2002al) ' selected a complete sample of HLIRGs in a 
manner independent of obscuration, inclination or AGN content 
and included sub-mm data (sub-mm data introduce a tight con- 
straint on the SB luminosities), finding that all HLIRGs in the 
sample were composite objects. 



In this paper we present a study of HLIRGs SED with two 
majors improvements and one limitation compared with the ear- 
lier studies commented above: (a) we have greatly enlarged the 
wavelength coverage, from radio to X-rays, and (b) we have sig- 
nificantly increased the photometric data coverage. However, as 
a self consistent analytical model able to reproduce the whole 
SED at so broad frequency range would be very complex to com- 
pute (and beyond the scope of this paper), we have compared 
our constructed SEDs with empirical AGN and SB templates, 
instead of using analytical RTM as in previous studies. 

The paper is organized as follows. Section|2]briefly describes 
the HLIRG sample. Section [3] explains how we built the SED 
and the data used to this end, and Sect. |4] the methods we have 
employed to model the SEDs. Results are presented in Sect. |5] 
compared with previous studies of HLIRGs in Sect. |6] and dis- 
cussed in Sect.|7] Section[8]summarizes our conclusions. 

The Wilkinson Microwave Anisotropy Probe ( WMAP) con- 
cordance cosmology has been adopted througho ut this paper: 
Ho ^ 70 km s-i Mpc"', Q„, = 0.27, Qa = 0.73 (ISpergel et al.l 
I2OOI . 

2. The HLIRG sample 

The sa mple studi ed here is the one investigated in iRuiz et al] 
(12007). From the iRowan-RobinsonI (l2000l) sample of HLIRGs 
we selected those sources with public data available in the 
XMM-Newton Science Archive as of December 2004, and we 
added our own XMM-Newton AO-5 observations. 

We limited this sample to sources with redshift less than 
~ 2 to avoid strong biasing towards high redshift quasars. 
Nevertheless, selecting the sample by using the availability of 
X-ray data probably introduces a selection effect in favour of the 
presence of an AGN. We also rejected one source from the orig- 
inal sample, IRAS 13279H-3401. Using recent optical and MIR 
spectra, we have determined its redshift to be z ~ 0.02 (see 
Appendix A), much lower than the on e presented in the litera- 
ture (z - 0.36. lRowan-Robinsonll2000l) . Therefore, our estimate 
of its IR luminosity is 3 x 10'" Lq, even below that necessary 
to classify it as a LIRG. Hence, we have thirteen objects in our 
final sample (see Table[T]i. 

According to their optical spectra (derived from the litera- 
ture), two sources are classified as starburst galaxies and twelve 
sources present AGN features. Among the latter eight are classi- 
fied as 'type I', and four of them as 'type IT. All type II and 
one NL-SB gala xy are Compton-Thick (CT) candidates. See 
iRuiz et al ] (120071) for a further discussion on this sample. 

3. Data compilation 

Our goal is to construct a well sampled SED for each object in 
a broad frequency range, from radio to X-rays. To this end, we 
have carefully searched in the literature and in several astronom- 
ical databases. See AppendixiBlfor a complete description on the 
origin of the photometry data for each HLIRG. 

All data included in the SEDs (presented in Tables B.1-B.13, 
see Appendix |B]i have been converted to monochromatic flux 
density units, corrected for the Galactic reddening and blue- 
shifted to rest-frame. 



3.1. Radio 

Most of the HLIRGs in the sample have at least one observation 
in the radio range. These data come from different observations 
by VLA, ATCA, IRAM and other radio-telescopes. 
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3.2. Infrared 

Our sources are well observed in the IR band. There are 

photometric data from IRAS (Po int S ource Catalogue, 

[joint IRAS Science Working Grou3 119881: Faint Source 
Catalogue. iMoshir et al.lll99 0) or ISO for all the objects. Most 
of the m also have been o bserved with SCUBA in the sub-mm 
band dFarrah et al.ll2002ah . 

In addition, there are public Spitzer MIR data for several 
sources: IRAC photometric data and IRS spectra. We have re- 
duced the IRAC data and we made our own photometric mea- 
surements. We have re-binned the IRS spectra of our HLIRG in 
broad bands, avoiding known emission and absorption features 
(a further analysis of these MIR spectra will be presented in Ruiz 
et al., in preparation). Mo st of these sources also have NIR data 
from the 2MASS surve>Q (ICutri et al.ll2003h . 

3.3. Optical and UV 

Most of the optical data were obtained from the Sloan 
Digital Sky Survey-Data Release 50 (SDSS-DR5, 
lAdelman-Mc Carthv et all 12007 ) and SuperCOSMOS Sky 
SurvejQ (SSS). A few data in V and B bands were taken from 
the XMM-Newton Optical Monitor (OM). 

We have only a few data in the UV range, mostly from the 
OM. Other data come from lUE and FUSE observations. 

3.4. X-ray 

The X MM-Newton spectra previously studied in iRuiz et al] 
(l2007h are available. We have correct ed each X-ray spectrum for 
the line of sight Galactic absorption (iDickev & Lockmanll 1 990h 
and we have re-binned the data in just a few energy bandfl In 
addition, the X-ray and the OM data come from simultaneous 
observations, allowing us to check any variability effects. 

3.5. Overall description of the SEDs 

Figure [U shows the SEDs we have built for our sources. We have 
divided the sample in two classes accordingly to their optical 
spectral classification. On one hand we grouped objects classi- 
fied as type I AGN (named class A sources) and on the other 
hand objects classified as type II AGN and SB (named class B 
sources). 

From a purely phenomenological point of view, class A and 
B sources seem to show a different SED shape. Class A objects 
have a SED approximately flat from the FIR to the optical range 
(the typical shape of quasars' SED), while class B objects show 
a prominent broad IR bump dominating the emission over the 
rest of the spectrum. 

To check if the above distinction holds quantitatively, we 
compared the distribution of X-ray-to-IR and optical-to-IR flux 
ratios for class A and class B sources. We estimated the 
monochromatic fluxes at three different rest-frame wavelengths, 
in the IR (30 yum), optical (4400 A) and X-rays (2 keV) through a 
linear interpolation of the SED (these points lie in well-sampled 

' http://www.ipac.caltech.edu/2mass/ 

^ http://www.sdss.org/dr5 

^ http://www-wfau.roe.ac.uk/sss/ 

* Through our X-ray data reduction we did not detect the source 
IRAS 14026+4341. Even so, this sou rce has a counterpart in the 
2XMMi catalogue jWatson et alj|2009l) . We have considered the five 
energy band fluxes as in the 2XMMi catalogue. 



regions of the SEDs, so these are reasonable estimates of the 
continua at those energies). Fig |2] shows the distribution of the 
X-ray-to-IR {Fx/ Fir) and optical-to-IR (Fopt/^^iR) flux ratios 
for the class A (blue histogram) and class B (pink histogram) 
sources. The distributions seem to be different for both classes 
of HLIRGs. By using a Kolmogorov-Smirnov test, the probabil- 
ity that class A and class B samples come from different parent 
populations is 92.6% for the Fx/Fir distribution and ~ 99.7% 
for the Fopt/FiR distribution. 

This rough analysis of the SED properties is clearly limited, 
but the results seem to support our classification of HLIRGs in 
two classes. We suggest that, since the SED classification is di- 
rectly related to the optical spectra classification, the distinct 
SED shape of HLIRGs could be explained by different levels 
of obscuration in the line of sight and/or the relative contribution 
of the SB emission to the total output. 

4. SED fitting 

Once all the SEDs were built, our aim was to check for the 
presence of AGN and/or SB emission in these sources and es- 
timate the contribution of these components to the total output. 
We fitted all SEDs by using the minimization technique with 
a simple model based on the use of convenient templates (see 
Sect. 14. II for details). The fitting procedure and the SED tem- 
plates were i mplemented usi ng the modelling and fitting tool 
Sherpa (Free man et al.ll200l1) . included in the software package 
CI AO 3.40. 

Our model comprises two additive components, one associ- 
ated to the AGN emission and the other associated to the SB 
emission. We can express this model as follows: 

i^v = ^'BOL(«M^™+(l-«)Mf), (1) 

where Fbol is the total bolometric flux, a is the relative contri- 
bution of the AGN to F^ol, Fy is the total flux at the frequency v, 
while u^'^^ and uf^ are the normalized AGN and SB templates 
(i.e., the value of the integral over the whole range of frequencies 
is unity for each SED template). This model contains only two 
free parameters, F^ol (the normalization) and a. The bolometric 
luminosity can be estimated as Lbol - ^-ttD^Fsol, where Dl is 
the luminosity distance. 

The model we are adopting to fit the SED is somehow rough 
and does not provide a precise description of the SED features, 
so we expect a poor fit in terms of value. However, the en- 
tire SED shape, from the radio to soft gamma rays, depends on 
a large number of physical parameters which produce different 
SED shapes even among the same class of sources (AGN, SB, 
etc.). Moreover, the impact of the different individual physical 
quantities on the overall SED and, perhaps most importantly, 
the effect of their interplay and interaction on the overall SED 
shape is far from being robustly settled from a theoretical point 
of view. The development of an analytical or semi-analytical 
model would be of great importance, but given that such mod- 
els are difficult to build and likely not unique, they clearly are 
beyond the scope of this work. We propose instead the simpler 
template-fitting approach to discriminate, as a zeroth-order ap- 
proximation, the relative component contribution (AGN and/or 
SB) to the overall bolometric luminosity of each source. 

We have chosen the fit with the lowest reduced as our 
"best fit" model. As we said above, the value of ;i'^/d.o.f. >> 1 
even for these best fits. Nevertheless, this quantity varies sig- 
nificantly for most sources between the different combinations 

' http://cxc.harvard.edu/ciao3. 4/ 
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of templates we tested during the minimization process. In 
those objects where different types of templates obtained simi- 
lar values, we have chosen the template most consistent with 
previous results in the literature. 

Our templates were chosen to minimize the contribution 
of the host galaxy's non-SB stellar emission (see Sect. 14.11 ). 
but there could still be a remnant of this emission in the tem- 
plates. Therefore, by adding two different templates we could 
have summed twice this stellar emission. We checked this effect 
adding a stellar template to the mode0. The normalization of this 
component was free and negative, in order to subtract the "sec- 
ond" stellar contribution. The addition of the new component did 
not change the final results of the SED fitting, so we can reject 
any important stellar contamination in our templates. 

4.1. Templates 

The templates we have employed to model the SEDs of our 
sources are empirical SEDs of well observed SB and Seyfert 
(Sy) galaxies in the local universe (see Table|2]l. 

To reproduce the AGN contribution we used six AGN tem- 
plates: 



1 . Two mean SED of radio quiet local QSO (Fig. |3(a)|): a lu- 
minosity independent SED (.Elvis et al. ,1994; Richards et al.l 
[2006) and a luminosity-dependent one (IHop kins et al. 2007). 
The latter template is similar to the standard SED of QSO 
from Elvis et al. (1994), but the value of apx depends on 
the bolometric luminosity dSteffen et al.ll2006h . and the X- 
ray emission above 0.5 keV is modelled by a power law (F - 
1 .8) with a cut-off at 500 keV and a reflection component 
generated with the PEXRAV model (Magd ziarz & Zdziarskil 
Q995). Therefore, this template has two parameters: normal- 
ization (the bolometric flux of the AGN) and redshift. For a 
given flux and redshift, the bolometric luminosity is calcu- 
lated and, hence, the value of aox- The first parameter was 
left free to vary during the fitting, while the second was fixed 
accordingly to the redshifts obtained in the literature. 

2. Four Sy2 galaxies (Fig. |3(b)[ ): these objects have hydrogen 
column densities (A^h) varying from 10^^ cm"^ (Compton 
thin objects) to greater than 10^^ cm~^ (Compton thick ob- 
jects). They were selected from a sample of Sy2 galaxie s 
with minimal starburst contribution ( Bianchi et al.l |20()6[) . 
The AGN templates show two bumps, in the FIR and in 
the NIR-optical, except for the AGN3 template, which only 
present a broad IR bump. The differences between them are 
the relative height of these bumps, the position of their peaks 
and the ratio between the optical and X-ray fluxes. 

To represent the SB contribution we have chosen a set of 
four starburst galaxies well observed in the full spectral range 
(Fig. |4(a)[ ). We have tried to cover a broad range of burst ages, 
dust contents and SFR. These physical properties are reflected 
in the SEDs showing different levels of obscuration, width and 
wavelength peaks. 

1 . NGC 5253 is a low-metallicity star-forming dwarf galaxy. Its 
nucleus is the site of a highly obscu red and extremely young 
(< 10 Myr) burst of star formation dBeck et alJI 19961) . with a 
SFR~ 8 Mq yr-i . 

2. NGC 7714 is a young unobscured SB (Brandl et al."2004) 
with SFR~ 6 M(T) yr~' and a burst age between 3-5 Myr 
dGonzalez-Deleado e"tanil995h . 

* The SED of the elliptical galaxy M 87 was employed to model the 
stellar emission. 



M82 is an evolved pur e SB galaxy with SFR~ 10 Mq yr ' 
(Stric kland e"tani2004 . 

IRAS 12112H-0305 is a bright ULIRG powe red by SB and 
with severe Umits to any AGN contribution dimanishi et alJ 
120071: iNardini et al.ll2008l). The estimated SF R for this object 
is ~ 600 vr"' dFranceschini et al.ll2003l) . 



All of them show two bumps, peaking in the FIR and in the 
NIR-optical. The main difference between the templates is the 
relative height between these bumps and their widths. 

We included four SED templates built from sources which 
they harbour both an AGN and a SB (Fig. |4(b)| l: 

1. NGC 1068 is a Sy2 galaxy with a composite nature , i.e. it 
harbours a heavily buried AGN (A^h > 10^ cmj^. Mat t et al.l 
[1997!) and also an intense SB (iTelesco et al.fl984) . The bolo- 
metric luminosity of this object is roughly evenly divided 
between the two component. The SB emission dominates 
longward of 30 fim and the AGN dominates shortward of 
20-10;um. 

2. Mrk 231 is an ULIRG (Lir = 3.2 x 10'^) optica lly classified 
as a Broad Absorption Line QSO (lBertall2005l) with a mas- 
sive young nuclear SB which is res ponsible for 25%-4 0% of 
the nuclear bolometric luminosity (iDavies et al.ll2004l) . 

3. IRAS 19254-7245, the "Superantennae", is a double- 
nucleated ULIRG optically classified as a Sy2 galaxy, with 
intense star formatio n. The AGN cont ribution to the total 
output is ~ 40 - 50% dBerta et al.ll2003l). 

4. IRAS 22491-1808 is a Sy2 ULIRG (Berta 2005) where the 
AGN emission is ~ 70% of the bolometric luminosity 
dFarrah et al.ll2003l) . 

We fitted these composite templates to those HLIRGs where 
the initial AGN+SB model was insufficient to reproduce the data 
(see Sect. 15.2b . 

We extracted the photometric data for the templates us- 
ing VOSECQ and VOSpecfl sof tware. These utilities use Virtual 
Observatory dOuinn et al.l2004t) tools to extract photometric and 
spectral data from several astronomical archives and catalogues. 
The templates were improved with data from NED database 
in wavelength ranges where VOSED and VOSpec provided no 
data. These objects are well observed at all the frequency ranges, 
particularly in the NIR and optical bands. We rejected some re- 
dundant data and we tried to extract only the nuclear emission 
to avoid as much contamination from the host galaxy as possi- 
ble. To this end we have chosen only those data with a roughly 
constant aperture within the nucleus of the galaxy. 



5. Results 

Figures|5]and|6]show the SEE0and the best fit model selected for 
each object, and Table [1] summarizes the results of our analysis. 
See Sect. l5.4l for comments on some particular sources. 



5.1. Class A HLIRGs 

We have shown that our simple two-component SED model is 
a fair approximation for most of these HLIRGs (see Fig. |l(a)| i. 
We found that all class A HLIRGs but one (IRAS 14026-H4341, 



' http://sdc.laefF.inta.es/vosed 
http://esavo.esa.int/vospec 

' Several photometric points are upper limits. The most conservative 
approach was chosen for the fit. We set the point to zero and the upper 
error bar to the upper limit value. 
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see Sect. 15 .41 below) are well fitted with type I AGN templates, 
consistent with their optical classification, and an additional SB 
component is required in three objects. The AGN component 
dominates the bolometric output in four out of these six sources, 
while two objects present a powerful SB component, with 60%- 
70% contribution of the bolometric luminosity. 

We have then 3 sources with a SB-dominated 
SED (IRAS F12509+3 122, IRAS 14026 -h4341 and 
IRASF14218+3845), one AGN-dominated source with an 
important SB contribution (IRAS 18216-1-6418), and 3 objects 
which seem to be extremely luminous quasars with no particular 
differences fr om the local one s, judging from their SEDs and 
X-ray spectra jRuiz et al.ll2007t) . 

A noticeable result for the class A HLIRGs is that the AGNl 
template over-predicts the X-ray flux of these sources, as found 
in our previous X-ray analysis. These discrepancies in the X- 
ray band can not be related to variability effects, since the OM 
data, simultaneous to the X-ray observations, match well with 
other optical and UV data obtained in different epochs. When we 
modelled these objects with the luminosity-dependent AGNl-L 
SED template, we found a significant improvement in the fit in 
terms of for most sources (4 out of 6) and the X-ray emission 
is better predicted. This result is consistent with the known am 
luminosity rela tionship (Strateva et al. 2005; Steffen et al. 200q 
iKellv et al.ll20b8h . 

We must also note that the IR-to-bolometric ratio of these 
sources is within ~ 40 - 70%, so an important fraction of their 
bolometric output is not emitted in the IR range. Hence, strictly 
speaking, they should not be considered as HLIRGs, particu- 
larly those with a completely AGN-dominated SED, where less 
than 50% of their bolometric luminosity is in the IR. This "con- 
tamination" can b e expected given the selection criteria of the 
iRowan-RobinsonI (2000) parent sample, which simply selected 
those known sources with Lir > 10'^ Lq. 

5.2. C\ass B HLIRGs 

We found that these sources are fitted with a dominant SB com- 
ponent and, in most cases, a minor AGN contribution (< 10%). 
However our model present some problems for class B HLIRGs 
that we did not find in class A objects (see Fig.|6l): 

1. The level of obscuration in the observed X-ray spectra is 
higher than the one expected from the AGN templates. 

2. Most sources show an excess in the MIR-NIR band not mod- 
elled by these templates, i.e. the width of the IR bump seems 
to be broader than the bumps in the starburst templates. 

3. The peak of the template does not match the IR peak of the 
SEDs in several sources. 

In order to improve the fit quality for the class B sources, we 
repeated the SED fitting using a set of templates from compos- 
ite sources (see Sect. 14.1b . where both AGN and SB emission 
are significant. By using these composite templates, we found 
that the statistical quality of our fits was significantly improved 
for all but one case (IRAS F15307+3252, see Sect. ISHbelow). 
For most objects, the obtained with any of the composite 
templates is significantly lower than the obtained with any 
combination of pure AGN and pure SB templates. CPl is the 
best fit template for 4 out of 6 sources, consistent with their 
spectral classification (type 2 AGN) and X-ray obscuration level 
(Compton-Thick). Two sources are best fitted with the CP2 tem- 
plate. 



IRAS F00235H-1024 is the only source that still shows a sig- 
nificant IR excess, which suggest that the SB contribution may 
be larger in this source that in the CPl template (~ 50%). 

5.3. Fitting without X-ray data 

In order to check how much X-ray data influence the SED fit- 
ting results, we exclude X-ray data from the SED fitting pro- 
cedure. Class A sources are still well represented by the same 
models (see Table [T] columns labeled as "no X-rays"), while 
class B galaxies are preferentially fitted with an AGN3 template 
(Compton thin model) and a SB component. Moreover, the AGN 
contribution grows significantly in most sources, particularly in 
the class B sources. When X-rays are included, a severe limit is 
imposed and the AGN contribution decreases dramatically. This 
shows that X-rays are important to obtain an accurate model with 
our technique and, hence, a better estimation of the contribution 
of each component to the total output. 

5.4. Notes on particular sources 
IRAS 14026-1-4341 

This source is optic ally classified as a type I AGN 
(iRowan-Robin son 2000), in agreement with the SDSS classifi- 
cation, and recent MIR Spitzer data also suggest the presence 
of an AGN in this object (Ruiz et al., in preparation), but our 
best fit model is obtained by using the SB2 template. The X-ray 
data impose a severe constraint, rejecting the AGN templates 
that predict a higher emission in the X-ray band. If we fit again 
this SED using no X-ray data (see Sect. 15. 3b we find that the best 
fit is obtained by AGN1-HSB2. 

The X-ray emission of this source seems to be affected by 
absorption (see Fig. |5(d)| l: it is not detected in the soft X-ray 
band (0.5-2 keV) and its 2XMMi hai'dness ratio {HR3 ~ - 0.2a 
is co nsistent with an X-ray absorbed AGN (Delia Cec a et al.l 
120041) . This indicates IRAS 14026-1-4341 as an X-ray absorbed 
QSO. These o bjects are of ten embedded in ultraluminous star- 
burst galaxies jPage et al] 12007). and they have been pointed 
out as a transitional phase in an evolutionary sequence relating 
the growth of massive black holes to the formation of galaxies 
dStevens et alJl2005HPage et al.ll2007[) . 

Under these circumstances, we have selected as best fit the 
model resulting from fitting the SED without X-ray data. We 
must note, however, that both models (pure SB or AGN-hSB) 
poorly fit the data between ~ 1-100 fim. The observed IR ex- 
cess, may be related to the X-ray emission absorbed and repro- 
cessed in the IR, can not be reproduced by AGNl template (an 
unabsorbed template). 

IRAS F1 5307-1-3252 

This object has been optically classified as a QSO 2 
(iRowan-RobinsonI |2000|) and there is strong evidence in X- 
rays favouring the presence of a heavily obscured AGN 
(llwasawa et alJ 2005h . However we have found that its SED best 
fit, in terms of x'^, is obtained with a SB template with minor 
AGNl contribution. The CPl template is also a fair fit, but with 
a slightly worse x^- 

Previous analyses of the IR emission of this HLIRG 
(iDeane & Trendiamll200ll IVerma et al1l2002h suggest that the 



= ^t^l^^^Z^^ ' -here CR is the count rate in the 
given energy band. 
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SB contribution is considerably lower than what we found us- 
ing a pure SB template. Hence, we have selected the CPl as 
"best fit", which is also consistent with its optical classification, 
to estimate the AGN and SB contribution to the bolometric lu- 
minosity. 

6. Comparison with previous results 

6.1. X-ray emission 

We can estimate the expected X-ray luminosity of the AGN and 
SB components for each source in our sample using the param- 
eters obtained in our SED analysis, and compare with the X-ray 
luminosities calculated through XMM-Newton observations. 

We have seen that the AGN SED of these sources is bet- 
ter modelled with a luminosity-dependent template. Hence, we 
have employed the relation obtained by Sani et al. (private com- 
municationo to estimate the intrinsic 2-10 keV luminosity for 
a given AGN bolometric luminosity: 

I- /T \ -0-357 

^±1^^0.043^] (2) 

LbOL \ 1045 ; 

Figure|7]shows those sources detected in X-rays and with an 
AGN component in their SED model. We plotted the bolometric 
luminosity of the AGN component versus the intrinsic (absorp- 
tion corrected ) 2-10 keV luminosity (see Table O, as calculated 
in lRuizet al.l(l2007l) . 

Most sources are scattered roughly following the Eq.|2]esti- 
mate. This scatter is probably related to the intrinsic dispersion 
in X-ray luminosities of AGN, i.e. for a given bolometric lu- 
minosity, there is a broad range of possible X-ray luminosities 
(ISteffenetalJl2006t) . 

There are, however, three sources (PG 1206-1-459, 
IRAS F12509-H3 122 and IRAS 14026+4341) with X-ray 
luminosities much lower than the estimated by Eq. |2l The 
X-ray luminosity of IRAS 14026+4341 was calculated using 
the 2XMMi X-ray fluxes so it is not corrected by absorption. 
Hence, this large discrepancy between the prediction and the 
observed luminosity is likely another sign of X-ray absorption 
(see Sect.l5H). 

For the other two sources iRuiz et alJ (l2007h did not find any 
sign of X-ray absorption. This effect could be, in principle, due 
to an overestimate of the AGN contribution to the bolometric 
luminosity. If we assume that the difference between the bolo- 
metric luminosity calculated using the SED fitting and that es- 
timated using Eq. |2]is completely originated due to star forma- 
tion, we find that the SB contribution to the total output should 
be larger than 90% in these two sources. Such a powerful SB 
must be clearly reflected in the SED shape, but we did not find 
this kind of deviation in the SED analysis of these sources. The 
X-ray weakness of these HLIRGs can not therefore be related to 
the underestimate of the SB contribution to the bolometric out- 
put, or due to X-ray absor ption. T hey seem to be intrinsically 
weak X-ray sources (Leigh lv et alj|20 01, 2007). 

6.2. IR SED: comparison witli previous work 

The IR (1-1000/im) SED of our sources has been previ- 
ously studied: Rowan-Robinson (2000), Farrahetal. (2002a) 

' ' This ratio is obtained from the lSteffen et al.l l l2006h relation between 
X-ray and 2500 A luminosities and then linking the 2 500 A luminosity 
with the bolometric one through the lElvis et alJ ( ll994h SED. 



and IVerma et al.l (l2002h modelled it using RTM. We estimated 
the IR luminosities of our models, integrating between 1- 
1000 jum, and compared their results with ours (see Table[3]l. 

The IR luminosities estimated through our SED fitting and 
that estimated using RTM match fairly well (see Fig. |8(a)[ ) 
for most sources. For three objects, our luminosity estima- 
tion is almost an order of magnitude greater than the RTM 
estimation, probably because our best-fit models overestimate 
the FIR-submm emission (see Figs. |5(b)[ |5(f)| and |6(f)[ ). This 
spectral emission is problably better recovered by using RTM. 
Nevertheless, in spite of this large disagreement in luminosities, 
our AGN contribution estimates are consistent with those ob- 
tained through RTM, as Fig. [8(b)] shows. 

The latter plot shows that our AGN contribution estimates for 
most sources are roughly consistent with those obtained through 
RTM. We can conclude that our simple model based on tem- 
plates is a fair method to obtain a first estimate of the AGN and 
SB relative contribution to the IR output. 

7. Discussion 

The broad band SEDs of the HLIRGs presented in this work 
can be roughly well fitted using templates, and their best fits 
are consistent with the optical classification of most sources (9 
out of 13). Among class A sources we found three objects fitted 
with pure type 1 AGN templates. They seem to be very luminous 
quasars and, since most of their bolometric output is not emit- 
ted in the IR band, should not be considered as proper HLIRGs. 
Four out of seven class A HLIRGs require, in addition to a type 
1 AGN template, a SB component which is, in three cases, domi- 
nant with respect to the AGN. The AGN emission in four sources 
is consistent with a luminosity-dependent SED. 

On the other hand, we have found that class B sources can 
not be fitted with simple combination of pure AGN and pure SB 
templates: a composite template is needed, where AGN and SB 
phenomena are both significant. This suggests that there should 
be some feedback between accretion process and star formation 
that changes the shape of the SED in a way that can not be imi- 
tated just by combining a pure SB and a pure AGN components. 
The main observational imprint of this feedback seems to be an 
excess in the SED around ~ 10 yum with respect to the predicted 
emission of pure AGN and pure SB combined model. 

Our division between class A and class B sources is based 
on the optical spectral classification, and since all objects show 
a significant AGN emission, it seems that the SED shape dif- 
ferences between the two groups c ould be an i nclination effect 
as in the unified model of AGNs dAntonucci & Miller 1985|: 
those HLIRGs where we have a direct view of the nucleus are 
luminous QSO and show a class A SED, while those HLIRGs 
seen through the molecular torus and/or other obscuring mate- 
rial show a class B SED. The comparable mean SB contribution 
of class A (excluding pure AGN sources) and class B sources is 
consistent with this hypothesis. Within this scenario, all types of 
HLIRGs belong to the same class of sources, seen at different 
inc hnation angle s. 

iFarrah etaP (l2002al) proposed, however, that HLIRG popu- 
lation is comprised of (1) mergers between gas rich galaxies, as 
found in the ULIRG population, and (2) young active galaxies 
going through their maximal star formation periods whilst har- 
bouring an AGN. 

The A^H distribution we found in the X-ray study seems 
to favour the two population hypothesis. In a pure inclination 
scenario we would expect a broad range of X-ray absorption, 
from no absorbed to heavily absorbed sources. However we 
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found only objects with no significant intrinsic absorption (all 
but one class A sources) or CT absorbed objects (all class B 
sources). Since AGNs observed in ULIR Gs usi ially show heavy 
absorption in X- rays dFranceschini et al.ll2003l:IPtak et al.ll2003t 
iTenp et al.l2005h . in principle class B sources could represent the 
high luminosity tail of ULIRG population, while the strong SB 
found in class A HLIRGs could represent young active galaxies 
experiencing their maximal star formation, without being in in- 
teracting systems (i.e. with little connection with a recent major 
merger). 

The study of the host galaxy morphology and environ- 
ment of HL IRGs also support the two population hypothesis. 
iFarrah et al.l ( 12002b ) found in a sample of nine HLIRGs ob- 
served by HST both strongly interactive systems and objects 
with no clear signs of ongoing interactions. Five sources of 
this sample are also included in ours: IRAS F00235+ 1024 and 
IRAS F15307 -1-3252 (class B objects) show signs of strong in- 
teractions, while IRAS F12509-H3 122, IRAS F14218-f3845 and 
IRAS 16347-1-7037 (class A objects) are isolated systems. This 
result favours our suggestion that class B HLIRGs could be ob- 
jects in the extreme bright end of the ULIRG population distri- 
bution. 

Hence, while class B HLIRGs share common properties with 
ULIRGs (high levels of X-ray obscuration, strong star forma- 
tion, signs of mergers and interactions), class A HLIRGs seem 
to be a different class of objects. Excluding the 3 pure AGN 
sources, class A objects could be am ong the young active galax- 
ies proposed bv lFarrah et al.l(l2002al) . The powerful SB we found 
in these sources, and the large amo unts of gas available to fuel 
the star formation (as calculated bv iFarrah et al.ll20 02a). along 
with the non detection of mergers or interactions in these sys- 
tems, support this idea. Moreover, the SB emission of the bona 
fide class A HLIRGs is modelled with young SB templates (SB 1 
and SB2) in all but one object (IRAS 18216-H6418), which is 
modelled with an old SB (SB3). This source could be a more 
evolved object. 

Therefore, sources in our sample likely belong to three dif- 
ferent populations: 

1 . Very luminous QSO with minor star formation activity. 

2. Young, isolated active galaxies undergoing their first episode 
of major star formation with Uttle connection with a recent 
major merger. 

3. Galaxies which have recently experienced a 
merger/disturbance that brought lots of gas and dust in 
the inner regions. This event trigger both the star formation 
and the AGN activity in a heavily obscured environment. 
These objects suit well as the high luminosity tail of the 
ULIRG population. 

Nevertheless, our sample of HLIRGs is not complete in any 
sense and we cannot derive further conclusions about the global 
properties of the HLIRG population. Further studies established 
on larger and complete samples of HLIRGs are needed to con- 
clude if the division between class A and class B objects is just 
due to an inclination effect, or is based on intrinsic differences 
of their physical properties. 

8. Conclusions 

In this paper we have built and analysed the multi-wavelength 
SED (from radio to X-rays) of a sa mple of 13 HLIR Gs, pre- 
viously studied in detail in X-rays (iRuiz et al.ll2007h . We as- 
sembled the SEDs using public data in several astronomical 



databases and in the literature, and we modelled them using 
templates. Most sources are roughly well fitted with this simple 
model and we find AGN relative contributions consistent with 
those inferred by previous analyses of the IR SEDs of HLIRGs 
using radiative transfer models. 

We divided the HLIRGs in two groups, accordingly to their 
optical spectral classification: class A (type 1 AGNs) and class 
B (type 2 AGNs and SB) sources. A first look at their SED 
shape indicates some differences between the two classes: class 
A sources show a roughly flat SED between the IR and the opti- 
cal, while class B sources have a prominent IR bump dominating 
the rest of the emission. 

A significant fraction (3 out of 7) of class A HLIRGs seem to 
be very luminous quasars with no particular deviations from the 
local quasars. Strictly speaking these objects should not be con- 
sidered HLIRGs since most of their bolometric output is emitted 
outside the IR band. The SED of these QSO are consistent with 
a luminosity-dependent quasar template. The remaining class A 
sources show significant additional SB components, which are 
dominant in all but one object. Given their strong SB activity and 
the lack of any sign of mergers in these systems, these HLIRGs 
could be very young galaxies experiencing their first episode of 
maximal star formation. 

Class B HLIRGs show an IR excess that can not be mod- 
elled with any combination of our selected pure AGN and pure 
SB templates. This feature can be properly fitted using compos- 
ite templates (SEDs from objects where AGN and SB emission 
are both important). This suggests that a significant fraction of 
the emission of this class of objects is originated in a SB. This 
also shows that the feedback between accretion and star forma- 
tion processes modifies the SED of class B HLIRGs in a way that 
can not be replicated by just the addition of pure AGN and pure 
SB independent templates. Class B HLIRGs share many proper- 
ties with ULIRGs (high X-ray absorption, strong star formation, 
signs of mergers and interactions), so they could be just the high 
luminosity tail of this population. 

Therefore, we have found some evidence supporting the idea 
that bona fide HLIRGs are composed of two populations: young 
active galaxies with no sign for recent mergers most likely going 
through their first episode of strong star formation, and the high- 
luminosity end of the ULIRG population, where both the SB and 
AGN are likely triggered by a recent merger/interaction. Further 
observational studies based on larger and, most importantly, 
complete samples of HLIRGs are needed to obtain stronger ev- 
idence for this hypothesis. Moreover our simple template-fitting 
approach should be complemented with RTMs (or other theoret- 
ical models of AGN and SB emission), since the two approaches 
are complementary in many ways and their combination may 
shed further light onto the relative SB-AGN contribution and on 
the feedback processes that take place in the most interesting 
HLIRGs, namely those that are well represented by composite 
templates within our approach. 
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Table 1. Best fit models for the HLIRG's SEDs 
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Source 


z 


Type 


CT" 








Best Fit model* 






log Lbol ' 


AGN / SB"* 










all data'' 




no X-rays' 




composite temp.* 


(ergs"') 












Model 




a 


Model 


a 








Class A HLIRGs 


PG 1206+459 


1.158 


QSO 


N 


AGNl-L 




1 


AGNl 




1 


AGNl-L 


48.4 


1/0 


PG 1247+267 


2.038 


QSO 


N 


AGNl-L 




1 


AGNl 




1 


AGNl-L 


49.2 


1/0 


IRAS F12509+3 122 


0.780 


QSO 


N 


AGNl-L 


SBl 


0.3 


AGNl 


SB4 


0.5 


AGNI-L-hSBI 


47.7 


0.3/0.7 


IRAS 14026+4341 


0.323 


QSO 


N 




SB2 





AGNl 


SB2 


0.3 


SB2 


46.7 


0.3/0.7 


IRASF14218+3845 


1.21 


QSO 


N 


AGNl 


SBl 


0.4 


AGNl 


SBl 


0.3 


AGNl+SBl 


47.2 


0.4/0.6 


IRAS 16347+7037 


1.334 


QSO 


N 


AGNl-L 




1 


AGNl 




1 


AGNl-L 


48.9 


1/0 


IRAS 18216+6418 


0.297 


QSO 


N 


AGNl 


SB3 


0.8 


AGNl 


SB3 


0.8 


AGN1+SB3 


47.4 


0.8/0.2 


Class B HLIRGs 


IRAS F00235+ 1024 


0.575 


NL-SB 


Y 




SB3 







SB3 





CPl 


46.7 


-0.5 / -0.5 


IRAS 07380-2342 


0.292 


SB 


N 


AGN4 


SBl 


0.06 


AGN3 


SBl 


0.3 


CPl 


47.0 


-0.5 / -0.5 


IRAS 00182-7112 


0.327 


QSO 2 


Y 


AGN3 


SB4 


0.06 


AGN3 


SB3 


0.3 


CPl 


46.6 


-0.5 / -0.5 


IRAS 09104+4109 


0.442 


QSO 2 


Y 


AGN4 


SB4 


0.09 


AGN3 


SBl 


0.8 


CP2 


47.3 


-0.7 / -0.3 


IRAS 12514+1027 


0.32 


Sy2 


Y 


AGN5 


SB4 


0.06 


AGN3 


SB2 


0.9 


CP2 


46.7 


-0.7/ -0.3 


IRASF15307+3252 


0.926 


QSO 2 


Y 


AGNl 


SB3 


0.03 


AGN3 


SBl 


0.8 


CPl 


47.9 


-0.5/ -0.5 



" Compton Thick candidates. 

' The best fit adopted to estimate the bolometric luminosity and the AGN and SB fraction is marked in bold fonts. 

Bolometric luminosity in CGS units. 
'' Fraction of the bolometric luminosity originated in AGN and SB. Calculated through the parameter a of the best fit model. 
' Best fit using our original set of templates. 
^ Best fit not using X-ray data. 

* Best fit including the templates of composite sources. 
Table 2. SED templates used as models. 



Label 


Source 


Description 


AGNl 




local quasar's mean SED' 


AGNl-L 




luminosity-dependent QSO SED- 


AGN3 


NGC5506 


Sy2, A^H = 3 X 1022cm-2 


AGN4 


NGC4507 


Sy2, A^H = 4 X lO^'cm-^ 


AGN5 


Mnk3 


Sy2,AfH = 1.4 X lO^'^cm-^ 


AGN6 


NGC 3393 


Sy2, Nu>lx lO^^cm-^ 


SBl 


NGC 5253 


Young and dusty SB 


SB2 


NGC 7714 


Young and unobscured SB 


SB3 


M82 


Old SB 


SB4 


IRAS 12112+0305 


ULIRG 


CPl 


NGC 1068 


Composite template: AGN: ~ 50% 


CP2 


Mnk231 


Composite template: AGN: ~ 70% 


CP3 


IRAS 19254-7245 


Composite template: AGN: ~ 45% 


CP4 


IRAS 22491-1808 


Composite template: AGN: ~ 70% 



References. (D iRichards et al.ll2006l : (2) lHoDkins et alj2007l 
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Table 3. IR and X-ray luminosities. 



Source 




logLf™" 


logLfB" 


l02 L'°' * 


J ^AGN * 
^IR.RTM 


'"6 Mr.rtm 


logLf^'- 






(erg s-') 


(erg s-') 


(erg s-') 


(erg s-') 


(erg s-') 


(erg s-') 


(erg s-') 


(cm--) 


Class A HLlRGs 


PG 1206+459 


48.0 


48.0 





47.8 


47.8 


< 46.7 


45 11+""^ 

-QQ4 




PG 1247+267 


48.8 


48.8 





47.9 


47.9 


<46.8 


45.93+°-°^ 




IRAS F12509+3 122 


47.6 


46.8 


47.5 


47.0 


46.8 


46.6 


42 26+^-^^ 




IRAS 14026+4341'' 


46.5 


45.8 


46.5 


46.5 


46.3 


46.1 


42 7+f(i 




IRASF14218+3845 


47.1 


46.5 


46.9 


46.9 


46.1 


46.8 


44 60+'^-^3 




IRAS 16347+7037 


48.5 


48.5 





47.7 


47.7 


<46.8 


46 00+" '" 




IRAS 18216+6418 


47.1 


46.9 


46.6 


46.8 


46.6 


46.4 







Class B HLIRGs 



IRAS F00235+ 1024 


46.7 


46.4 


46.4 


46.7 


46.4 


46.4 


< 42.2 


> W 


IRAS 07380-2342 


47.0 


46.7 


46.7 


47.0 


46.8 


46.5 


< 41.7 




IRAS 00182-7112 


46.6 


46.3 


46.3 


46.7 


<46.5 


46.7 


44 82+""' 


> 10^5 


IRAS 09104+4109 


47.3 


47.1 


46.8 


46.8 


46.8 


< 46.2 


43. 3+'-'* 


> 10^5 


IRAS 12514+1027 


46.7 


46.5 


46.2 


46.5 


46.2 


46.2 


(4+f ) X 10^' 


IRASF15307+3252 


47.9 


47.6 


47.6 


46.9 


46.6 


46.7 


4"; AQ+l>'^9 


> 10^5 



a) IR luminosities (1-1000 pm) estimated using our SED fitting. 

b) IR luminosities (1-1000 yum) estimated by the analysis of the IR SED using RTM jRowan-Robinsonll 20001 : iFarrah et alll2002ah . 

c) Absorption corrected 2-10 keV luminosities from Ruiz et al. (2007). 

d) Intrinsic absorption estimated using X-ray spectra (Rui z et al..2007 ). 

e) The X-ray luminosity of this source has been calculated from 2XMMi fluxes iWatson et al.ll200^ . and it is not corrected from absorption. 




Fig. 1. Rest-frame spectral energy distributions of the sample. Fluxes are shifted for clarity. 
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1 




-3.0 -2.5 
X-ray-to-IR 

(a) Fx/F,R 




Fig. 2. Distribution of (a) X-ray-to-IR and (b) Optical-to-IR flux ratios for class A (dark grey, blue in the colour version) and class 
B (light grey, pink in the colour version) HLIRGs. 




10= 10' 10'° 10" 10>= 10'= 10'« 10'= 10'= 10" 10'= 10" 10= 10= 10'" 10" 10'= 10'= 10'* 10'= 10'= 10" 10'= 10'" 

V (Hz) V (Hz) 

(a) AGN I templates. (b) AGN II templates. 



Fig. 3 . AGN templates, (a) The top line (blue in flie colour version) is the standar d SEP for radio quie t quasar ( AGN 1, [Richards et al-l 
l2006l) . The group below is the luminosity-dependent SED for quasar (AGNl-L. lHopkins et al.ll2007h . plotted for several bolometric 
luminosities (the top line - red in the colour version - is for 10'" Lq and the bottom black line is for 10'^ Lg)- (b) Listed downwards: 
NGC5506 (AGN3,A?H = 3 x lO^^cm^^), NGC4507 (AGN4,A?h = 4x lO^^cm^^), Mrk3 (AGN5, A^h = 1.4x lO^'^cm-^), NGC 3393 
(AGN6, A^H > 1 X 10-^ cm^-). The SED fluxes are shifted for clarity. See Sect. ElUfor details. 




Fig. 4. (a) Pure starburst templates. Listed downwards: NGC5253, NGC7714, M82, IRAS 12112+0305. (b) Composite templates 
(AGN H- SB). Listed downwards: NGC 1068, Mrk231, IRAS 19254-7245, IRAS 22491-1808. The SED fluxes are shifted for clarity. 
See Sect.|4]T]for details. 
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Fig. 5. Rest-frame Spectral Energy Distributions of class A HLlRGs and their best fit models. Dotted lines (red in the colour version) 
are the AGN components and dashed lines (green in the colour version) are the SB components. Black solid lines are the sum of the 
AGN and SB components. 
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Fig. 6. Rest-frame Spectral Energy Distributions of class B HLIRGs and their best fit models. Symbols as in Fig. |5] The long-dashed 
lines (blue in the colour version) are the best fits obtained using composite templates (see Sects. l4~T] and l5.2l) . 
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Fig. 7. Bolometric versus observed, absorption-corrected 2-10 keV AGN luminosities. Squares (blue in the colour version) are class 
A HLIRGs, triangles (red in the colour version) are class B HLIRGs. The dotted line reflects the ratio between these luminosities 
obtained by Sani et al. 
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Fig. 8. (a) Total IR luminosity estimated using our templates (Lir) and using radiative transfer models (Lj^*^). (b) AGN to total IR 
luminosity ratios estimated through our model {Rkgh) ™d using radiative transfer models (^^gn-*' Symbols as in Fig.|7] The dotted 
lines mean equal values. 
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Appendix A: IRAS 13279+3401 

The object IRAS 13279+3401 has been previously classified 
as a QSO, and the IR luminosity estimated through the red- 
shift presented in the literature (z=0.36, Rowan-Robinson 2000; 
iFarrah et alj|2002al) exceed the HLIRG limit. However, we have 
now strong evidence showing that this source is a much closer 
galaxy. 

Figure A |l(a)| shows its optical spectrum obtained by the 
2.5m Isaac Newton Telescope, where we do not find any type 
I feature. A QSO with z=0.36 should present a broad emis- 



sion line at ~ 6600 A. We estimate z=0.023 for this spectrum, 
from stellar absorption features. 

We have also the MIR spectrum of this source, observed by 
Spitzer (see Fig. A |l(b)[ ). We hav e estimated the redsh ift of the 
source using a SB template from iNardini et al.] (l2008l ). We red- 
shifted the template matching the most important spectral fea- 
tures and we found z ~ 0.02, which is consistent with our esti- 
mate using the optical spectrum. The IR luminosity derived from 
this redshift is ~ 3 x 10'" Lq, well below the HLIRG Hmit and 
even below the LIRG. 
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4500 5000 5500 6000 6500 7000 7500 

Wavelength (A) ^ (M™) 

(a) Optical spectrum (2.5m INT). ^b) MIR spectrum {Spitzer Infrared Spectrograph). 

Fig. A.l. (a) Optical and (b) MIR spectra of IRAS 13279+3401 in the o bserver frame. The slashed hne (green in the colour version) 
in the right hand side plot is the SB template from lNardini et al.) (l2008l) . 
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Appendix B: Tables of data 

Along this appendix a table is presented for each object with the 
fluxes employed to build the SEDs and the origin of each data. 
The re-binned spectra from XMM-Newton and Spitzer are not 
included in these tables. Fluxes shown with no errors are upper 
limits. 

Table B.l, Photometric data for IRAS 00182-7112. 



v(Hz) 




Fy (Jy) 




Error (ly) 




Ref. 


8.43 X 


10« 


4.23 X 


10- 


-1 


1.28 X 


10- 


-2 


SUMSS (NED) 


1.40 X 


10' 


3.17 X 


10- 


-1 


3.00 X 


10- 




ATCA (Drake et alJ2004l) 


2.50 X 


10' 


1.97 X 


10- 


-1 


3.00 X 


10- 


-3 


ATCA ('Drake et al. 2004") 


4.80 X 


10' 


9.80 X 


10- 


-2 


3.00 X 


10- 


-3 


ATCA 1 Drake et al. 2004) 


8.60 X 


10' 


5.70 X 


10- 


-2 


3.00 X 


10- 


-3 


ATCA 1 Drake et alJ2004l) 


3.00 X 


10'2 


1.19 X 


10" 


1.19 X 


10- 


-1 


IRAS (NED) 


5.00 X 


10'2 


1.20 X 


10" 


8.37 X 


10- 


-2 


IRAS (NED) 


1.20 X 


10" 


1.33 X 


10- 


-1 


1.02 X 


10- 


-2 


IRAS (NED) 


2.50 X 


10'3 


6.02 X 


10- 


-2 








IRAS (NED) 


5.23 X 


10" 


2.23 X 


10- 


-2 


4.78 X 


10- 


-3 


Spitzer -IRAC 


8.44 X 


10" 


2.53 X 


10- 


-3 


7.47 X 


10- 


-4 


Spitzer -IRAC 


1.39 X 


10'^ 


7.22 X 


10- 


-4 


7.83 X 


10- 


-5 


2MASS 


1.80 X 


10'" 


4.37 X 


10- 


-4 


7.34 X 


10- 


-5 


2MASS 


2.43 X 


10"* 


2.68 X 


10- 


-4 


4.34 X 


10- 


-5 


2MASS 


3.33 X 


10"* 


2.02 X 


10- 


-4 


5.57 X 


10- 


-5 


SSS 


4.28 X 


10"* 


2.28 X 


10- 


-4 


6.30 X 


10- 


-5 


SSS 


6.81 X 


10"* 


2.97 X 


10- 


-5 


5.40 X 


10- 


-6 


XMM-WcwtoK -DM 


1.48 X 


10" 


1.90 X 


10- 


-5 


1.30 X 


10- 


-5 


XMM-A'cH'/o« -DM 



Table B.2. Photometric data for IRAS F00235+1024. 



v(Hz) 


Fy (Jy) 




Error (Jy) 




Ref. 


1.40 X 10' 


7 70 y 1 n- 


-3 


■s nn V 1 n- 


-4 


VI A ^■^TBT^^ 




A V in- 


-3 






oL.urSrt 1 ranan et ai/zuuzajj 


fi(S y in" 

U.UU A L\J 


9.67 X 10- 


-1 








1 66 X 1 n" 


8 n5 X in- 


-1 


3 93 X in- 


-1 




3.15 X 10" 


4.78 X 10- 


-1 


1.48 X 10- 


-1 


ISO (NED) 


5.00 X 10" 


4.28 X 10- 


-1 


5.56 X 10- 


-2 


IRAS (NED) 


2.00 X 10" 


6.75 X 10- 


-3 


2.14 X 10 


-3 


ISO (NED) 


4.44 X 10" 


9.20 X 10 


-4 


3.80 X 10- 


-4 


ISO (NED) 


1.39 X 10"* 


1.01 X 10- 


-4 


5.00 X 10- 


-6 


?Fan-ah et al.'2002al) 


3.33 X lO'"* 


8.00 X 10- 


-5 


1.00 X 10- 


-5 


( Fan-ah et al. 2002ai) 


4.28 X 10'" 


3.10X 10- 


-5 


7.00 X 10- 


-6 


USNO-BI.O 


6.81 X 10"* 


1.80 X 10- 


-5 


4.00 X 10- 


-6 


USNO-BI.O 


3.02 X 10" 


4.68 X 10- 


-10 






XMM-Newton -EPIC 


1.45 X 10'** 


1.04 X 10- 


-10 






XMM'Newton -EPIC 



Table B.3. Photometric data for IRAS 07380-2342. 



v(Hz) 


Fy (Jy) 


Error (Jy) 




Ref. 


3.53 X 10" 


2.68 X 10-2 


4.20 X 10- 


-3 


SCUBA IFarrah et alJ2002al) 


6.66 X 10" 


9.67 X 10^' 






SCUBA IFarrah et alJ2002aD 


3.00 X 10'^ 


3.55 X 10" 


2.84 X 10- 


-1 


IRAS (NED) 


5.00 X 10" 


1.17 X 10" 


9.36 X 10- 


-2 


IRAS (NED) 


1.20 X 10" 


8.00 X 10-' 


8.00 X 10- 


-2 


IRAS (NED) 


2.50 X 10'^^ 


4.84 X 10"' 


3.39 X 10 


-2 


IRAS (NED) 


1.39 X lO'-* 


4.08 X 10-3 


1.49 X 10- 


-4 


2MASS 


1.80 X 10'" 


1.35 X 10-3 


9.71 X 10- 


-5 


2MASS 


2.43 X lO'" 


4.74 X 10-" 


4.79 X 10- 


-5 


2MASS 


3.33 X 10'" 


1.62 X 10-" 


2.23 X 10- 


-5 


DENIS 


5.45 X 10'" 


3.80 X 10-"^ 


2.10 X 10- 


-5 


XMM-Newton -CM 


3.02 X 10" 


2.51 X 10-' 






XMM-Newton -EPIC 


1.45 X 10'^ 


1.87 X 10-' 






XMM-Newton -EPIC 
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Table B.4. Photometric data for IRAS 09104+4109. 



v(Hz) 




Fy (Jy) 




Error (ly) 




Ref. 


1.40 X 


10' 


6.88 X 


10- 


-3 


1.31 X 


10- 


-4 


V \l\ CiLJ ) 




in' 


6.00 X 


10- 


-3 


1.00 X 


10- 


-3 




4. on V 




1.80 X 


10 


-3 


3.00 X 


10- 


-4 


V LA CjU ) 


J.JJ A 




9.54 X 


10 


-3 








oL-UriA 1 ranan et di. zuuza) 


U.OU A 


in" 


7.28 X 


10- 


-2 








Ai^UrSA dranan et dijzuuzi^) 


fin V 
J .UU A 


ini2 


4.38 X 


10- 


-1 








TD A C ^'MtJr^'^ 


^ nn V 

J.UU A 


ini2 


5.25 X 


10- 


-1 


4.20 X 


10- 


-2 


TD A C ("MCT^"! 


1.20 X 


ini3 


3.34 X 


10- 


-1 


1.30 X 


10 


-2 


TD A Q /'MPT^'l 


1.26 X 


IQii 


3.33 X 


10- 


-1 


1.70 X 


10 


-2 


Vnit'^iir TD AC 


1.43 X 




2.70 X 


10- 


-1 


7.00 X 


10 


-2 


TDTtJ ik'lfinrmnn P«t nll1Q8RI\ 

iKir 1 ivicinrnann et aijiyooit 


2.50 X 


1013 


1.30 X 


10- 


-1 


3.11 X 


10 


-2 




2.97 X 


10" 


8.80 X 


10- 


-2 


1.70 X 


10 


-2 


TDTP 1 K"li=>intmnn f>t !il 1 1 QRRh 


5.23 X 


10" 


2.64 X 


10- 


-2 


7.11 X 


10 


-3 


^nitmr TD AC 


8.44 X 


10" 


4.74 X 


10- 


-3 


1.21 X 


10- 


-3 


-TR AC 


1.39 X 


10'" 


8.21 X 


10- 


-4 


8.61 X 


10- 


-5 




1.80 X 


lO'"* 


4.59 X 


10- 


-4 


6.24 X 


10 


-5 




2.43 X 


lO'"* 


2.97 X 


10- 


-4 


4.64 X 


10 


-5 




3.33 X 


10"' 


2.42 X 


10- 


-4 


6.69 X 


10- 


-5 


sss 


3.28 X 


10'" 


5.06 X 


10 


-4 


7.45 X 


10 


-6 


SDSS-DR5 


3.93 X 


10'" 


7.04 X 


10- 


-4 


3.42 X 


10 


-6 


SDSS-DR5 


4.28 X 


10'" 


2.12 X 


10- 


-4 


5.85 X 


10 


-5 


SSS 


4.81 X 


10'" 


2.11 X 


10- 


-4 


1.69 X 


10 


-6 


SDSS-DR5 


6.28 X 


10'" 


1.51 X 


10- 


-4 


1.26 X 


10 


-6 


SDSS-DR5 


6.81 X 


10'" 


1.08 X 


10- 


-4 


3.97 X 


10 


-5 


SSS 


8.47 X 


10'" 


1.01 X 


10 


-4 


3.16X 


10 


-6 


SDSS-DR5 


9.67 X 


10'" 


5.90 X 


10- 


-5 


6.00 X 


10 


-6 


XMM-Newton -CM 


1.48 X 


10'5 


7.50 X 


10- 


-5 


1.80 X 


10 


-5 


XMM-A'cwto« -CM 


7.86 X 


10"* 


1.06 X 


10- 


-7 


5.91 X 


10 


-8 


BeppoSAX (NED) 



Table B.5. Photometric data for PG 1206+459. 



v(Hz) 




Fy (ly) 




Error (ly) 




Ref. 


4.90 X 


10" 


1.20 X 


10- 


-4 








VLA (NED) 


2.30 X 


10" 


1.80 X 


10 


-3 


4.50 X 


10- 


-4 


IRAM (NED) 


1.72 X 


10'2 


1.89 X 


10 


-1 


3.78 X 


10- 


-2 


ISO (NED) 


2.93 X 


10'2 


3.53 X 


10 


-1 


7.06 X 


10- 


-2 


ISO (NED) 


4.93 X 


10'2 


2.57 X 


10 


-1 


5.14 X 


10- 


-2 


ISO (NED) 


1.20 X 


10'3 


1.13 X 


10 


-1 








IRAS (NED) 


1.43 X 


10'3 


6.40 X 


10 


-2 


1.28 X 


10- 


-2 


ISO (NED) 


2.34 X 


10'3 


2.30 X 


10 


-2 


4.60 X 


10- 


-3 


ISO (NED) 


2.50 X 


10'3 


2.07 X 


10 


-1 


3.60 X 


10 


-2 


IRAS (NED) 


6.17 X 


10'3 


3.00 X 


10 


-3 








ISO (NED) 


1.39 X 


10'" 


2.64 X 


10 


-3 


9.46 X 


10 


-5 


2MASS 


1.80 X 


10'" 


2.57 X 


10 


-3 


1.22 X 


10- 


-4 


2MASS 


2.43 X 


10'" 


2.43 X 


10 


-3 


7.76 X 


10- 


-5 


2MASS 


3.33 X 


10'" 


2.27 X 


10 


-3 


6.26 X 


10- 


-4 


SSS 


3.28 X 


10'" 


2.78 X 


10 


-3 


1.26 X 


10- 


-5 


SDSS-DR5 


3.93 X 


10'" 


2.94 X 


10 


-3 


9.76 X 


10- 


-6 


SDSS-DR5 


4.28 X 


10'" 


2.51 X 


10 


-3 


6.93 X 


10- 


-4 


SSS 


4.81 X 


10'" 


2.92 X 


10 


-3 


9.45 X 


10- 


-6 


SDSS-DR5 


5.45 X 


10'" 


2.62 X 


10 


-3 


3.00 X 


10- 


-5 


yMM'Newton -OM 


6.28 X 


10'" 


2.53 X 


10 


-3 


7.01 X 


10 


-6 


SDSS-DR5 


6.81 X 


10'" 


2.25 X 


10 


-3 


1.30 X 


10- 


-5 


yMM'Newton -OM 


8.47 X 


10'" 


2.20 X 


10 


-3 


9.06 X 


10- 


-6 


SDSS-DR5 


1.05 X 


10'5 


1.00 X 


10 


-3 


9.21 X 


10- 


-5 


lUE 


1.46 X 


10'5 


5.25 X 


10 


-4 


7.25 X 


10- 


-5 


lUE 


2.11 X 


10'5 


1.45 X 


10 


-4 


8.32 X 


10- 


-5 


lUE 



Table B.6. Photometric data for PG 1247+267. 



V (nz; 




(Jy) 




Error (ly) 




Ref. 


1.49 X 


10' 


1.17 X 


10- 










VLA (NED) 


4.90 X 


10' 


7.20 X 


10- 


-4 


8.00 X 


10- 




VLA (NED) 


1.49 X 


10'" 


1.51 X 


10- 




2.20 X 


10- 


-4 


VLA (NED) 


2.30 X 


10" 


2.10 X 


10- 










IRAM (NED) 


1.72 X 


10'^ 


1.50 X 


10- 


- 1 








ISO (NED) 


2.93 X 


10'^ 


1.74 X 


10- 


- 1 


3.48 X 


10- 


-2 


ISO (NED) 


4.93 X 


10'^ 


2.36 X 


10- 


-1 


4.72 X 


10- 


-2 


ISO (NED) 


1.20 X 


10'^ 


1.13 X 


10- 


- 1 








IRAS (NED) 


2.34 X 


10'^ 


3.00 X 


10- 


''^ 


6.00 X 


10- 




ISO (NED) 


3.91 X 


10" 


1.00 X 


10- 




2.00 X 


10- 




ISO (NED) 


6.17 X 


10" 


9.00 X 


10- 










ISO (NED) 


1.39 X 


10'" 


3.55 X 


10- 




1.36 X 


10 




2MASS 


1.80 X 


10'" 


2.99 X 


10- 




1.45 X 


10 




2MASS 


2.43 X 


10'" 


2.93 X 


10 




l.lOx 


10 




2MASS 


3.33 X 


10'" 


6.31 X 


10- 




1.74 X 


10 




SSS 


3.28 X 


10'" 


3.01 X 


10- 




1.37 X 


10- 




SDSS-DR5 


3.93 X 


10'" 


2.66 X 


10- 




9.69 X 


10 


-6 


SDSS-DR5 


4.28 X 


10'" 


4.03 X 


10- 


-3 


1.11 X 


10 


-3 


SSS 


4.81 X 


10'" 


2.26 X 


10- 


-3 


8.02 X 


10 


-6 


SDSS-DR5 


6.28 X 


10'" 


2.13 X 


10 


-3 


6.42 X 


10 


-6 


SDSS-DR5 


6.81 X 


10'" 


3.54 X 


10- 


-3 


1.30 X 


10 


-3 


SSS 


8.47 X 


10'" 


2.27 X 


10- 


-3 


9.92 X 


10 


-6 


SDSS-DR5 


9.67 X 


10'" 


1.17 X 


10- 


-3 


8.00 X 


10 


-6 


XMM-Newlon - 


1.25 X 


10'^ 


6.70 X 


10- 


-4 


2.00 X 


10 


-5 


XMM-Newton - 


1.48 X 


10'^ 


3.70 X 


10- 


-4 


2.00 X 


10 


-5 


XMM-Newton - 



Table B.7. Photometric data for IRAS F12509+3122. 



v(Hz) 




F. (Jy) 




Error (Jy) 




Ref. 


1.40 X 


10" 


1.76 X 


10 


-3 


1.25 X 


10 


-4 


VLA (NED) 


3.53 X 


10" 


9.23 X 


10 


-3 








SCUBA iFan-ah et al."2002a|) 


6.66 X 


10" 


3.33 X 


10 


-1 








SCUBA ( Fan-ah et aL2002i]) 


3.00 X 


10'^ 


6.75 X 


10 


-1 








IRAS (NED) 


5.00 X 


10'2 


2.18 X 


10 


-1 


4.36 X 


10- 


-2 


IRAS (NED) 


1.20 X 


10'3 


1.03 X 


10 


-1 


2.75 X 


10 


-2 


IRAS (NED) 


1.39 X 


10'" 


1.18 X 


10 


-3 


8.78 X 


10 


-5 


2MASS 


1.80 X 


10'" 


7.75 X 


10- 


-4 


8.07 X 


10- 


-5 


2MASS 


2.43 X 


10'" 


9.77 X 


10 


-4 


6.01 X 


10- 


-5 


2MASS 


3.33 X 


10'" 


8.60 X 


10 


-4 








( Fan-ah et aI.I2002iJ) 


3.28 X 


10'" 


8.73 X 


10 


-4 


7.57 X 


10- 


-6 


SDSS-DR5 


3.93 X 


10'" 


8.39 X 


10 


-4 


3.56 X 


10- 


-6 


SDSS-DR5 


4.28 X 


10'" 


7.20 X 


10 


-4 


7.00 X 


10 


-5 


APM 


4.81 X 


10'" 


8.39 X 


10 


-4 


3.12 X 


10 


-6 


SDSS-DR5 


5.45 X 


10'" 


6.60 X 


10 


-4 


1.60 X 


10- 


-5 


XMM-Newton -OM 


6.28 X 


10'" 


8.32 X 


10 


-4 


2.89 X 


10 


-6 


SDSS-DR5 


8.47 X 


10'" 


7.22 X 


10 


-4 


4.89 X 


10 


-6 


SDSS-DR5 


9.67 X 


10'" 


4.94 X 


10 


-4 


6.00 X 


10- 


-6 


XMM-Newton -OM 


1.25 X 


10'^ 


4.68 X 


10 


-4 


7.00 X 


10- 


-6 


XMM-Newton -OM 


1.33 X 


10'5 


7.83 X 


10 


-4 


9.94 X 


10 


-7 


HST - FOS 


1.48 X 


10'^ 


4.69 X 


10 


-4 


1.60 X 


10 


-5 


XMM-A'evvton -OM 


1.75 X 


10'5 


5.16X 


10- 


-4 


1.38 X 


10- 


-6 


HST - FOS 
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Table B.8. Photometric data for IRAS 12514+1027. 



Table B.IO. Photometric data for IRAS F14218+3845. 



v(Hz) 




Fr (Jy) 




Error (Jy) 




Ref. 






7.78 X 


10- 


-3 


1.46 X 


10 


-4 


VLA (JNbU) 


^ on V 




7.55 X 


10 


-1 


1.51 X 


10 


-1 




s no V 

J.UU X 


lU 


7.12 X 


10 


-1 


5.70 X 


10 


-2 




1 1C\ V 
i .ZU X 


lU 


1.90 X 


10 


-1 


1.58 X 


10 


-2 




9 so V 


lU 


6.32 X 


10 


-2 








TD A C ^'MC^^^ 


S 7'^ V 

J.Z,J A 


in" 


3.46 X 


10 


-2 


4.03 X 


10" 


-3 




O-H-H- A 




1.49 X 


10 


-2 


4.82 X 


10 


-3 




1 V 




2.57 X 


10 


-3 


1.26 X 


10 


-4 


ZlVIrtOO 


1 RC\ V 

i .OU A 


in'* 


8.45 X 


10 


-4 


8.57 X 


10 


-5 


ZlVIrtoo 


2.43 X 


10'-* 


4.07 X 


10 


-4 


5.40 X 


10 


-5 


2MASS 


3.33 X 


10"* 


2.88 X 


10 


-4 


7.96 X 


10 


-5 


sss 


3.28 X 


lo'-* 


5.66 X 


10 


-4 


1.15 X 


10 


-5 


SDSS-DR5 


3.93 X 


10"* 


3.39 X 


10 


-4 


2.96 X 


10 


-6 


SDSS-DR5 


4.28 X 


10'* 


2.46 X 


10 


-4 


6.81 X 


10 


-5 


SSS 


4.81 X 


10"* 


2.82 X 


10 


-4 


2.30 X 


10 


-6 


SDSS-DR5 


6.28 X 


10'* 


1.16X 


10 


-4 


1.50 X 


10 


-6 


SDSS-DR5 


6.81 X 


10"* 


5.78 X 


10 


-5 


2.13 X 


10 


-5 


SSS 


8.47 X 


10'* 


4.79 X 


10 


-5 


3.53 X 


10 


-6 


SDSS-DR5 



Table B.9. Photometi-ic data for IRAS 14026+4341. 



y (Hz) 




f V (Jy) 




Error (Jy) 




Ref. 




1.40 X 


10' 


1.59 X 


10- 


-3 


1.39 X 


10 


-4 


VLA (NED) 




3.53 X 


10" 


7.53 X 


10- 


-3 








SCUBA (Farrah et alJ2002al) 


6.66 X 


10" 


9.40 X 


10- 


-2 








SCUBA (Farrah et alJ2002al) 


3.00 X 


10'2 


9.94 X 


10- 


-1 


2.39 X 


10- 


-1 


IRAS (NED) 




5.00 X 


10'2 


6.22 X 


10- 


-1 


5.60 X 


10- 


-2 


IRAS (NED) 




1.20 X 


10" 


2.85 X 


10- 


-1 


1.41 X 


10- 


-2 


IRAS (NED) 




2.50 X 


10" 


1.18 X 


10- 


-1 


2.71 X 


10 


-2 


IRAS (NED) 




3.91 X 


10'3 


3.60 X 


10- 


-2 


4.00 X 


10 


-3 


ISO (NED) 




7.05 X 


10'3 


1.60 X 


10- 


-2 








ISO (NED) 




1.39 X 


10'* 


8.47 X 


10- 


-3 


2.29 X 


10- 


-4 


2MASS 




1.80 X 


10'* 


4.75 X 


10- 


-3 


1.57 X 


10 


-4 


2MASS 




2.43 X 


10'* 


3.33 X 


10- 


-3 


1.01 X 


10- 


-4 


2MASS 




3.33 X 


10'* 


2.26 X 


10- 


-3 


6.25 X 


10- 


-4 


SSS 




3.28 X 


10'* 


3.83 X 


10- 


-3 


2.05 X 


10- 


-5 


SDSS-DR5 




3.93 X 


10'* 


2.70 X 


10- 


-3 


1.23 X 


10- 


-5 


SDSS-DR5 




4.81 X 


10'* 


2.82 X 


10- 


-3 


1.17 X 


10- 


-5 


SDSS-DR5 




5.45 X 


10'* 


2.44 X 


10- 


-3 


1.90 X 


10- 


-5 


XMM'Newlon 


-CM 


6.28 X 


10'* 


2.09 X 


10 


-3 


9.05 X 


10 


-6 


SDSS-DR5 




8.47 X 


10'* 


1.30 X 


10- 


-3 


7.19 X 


10- 


-6 


SDSS-DR5 




9.67 X 


10'* 


7.70 X 


10- 


-4 


4.00 X 


10- 


-6 


XMM-Newton 


-DM 


1.25 X 


lO's 


2.65 X 


10 


-4 


7.00 X 


10- 


-6 


XMM-Newton 


-DM 


1.48 X 


10'5 


1.59 X 


10- 


-4 


9.00 X 


10- 


-6 


XMM-Newton 


-DM 


8.46 X 


10"> 


1.42 X 


10- 


-10 


8.39 X 


10- 


-10 


2XMMi 




1.81 X 


10" 


2.14 X 


10 


-10 


3.50 X 


10 


-10 


2XMMi 




3.63 X 


10" 


7.82 X 


10 


-10 


3.76 X 


10 


-10 


2XMMi 




7.86 X 


10" 


6.81 X 


10- 


-10 


3.85 X 


10- 


-10 


2XMMi 




1.99 X 


10"* 


5.49 X 


10 


-10 


5.17 X 


10 


-10 


2XMMi 





v(Hz) 


Fv (Jy) 




Error (Jy) 




Ref. 


j.Jj X lU 


O.DD X lU 


-3 






oL.UdA (rarran et aijzuuzai) 


O.OU A ixj 


A IV/ 


-1 






oL.uda (rarran et ai^zuuzai) 


J. I J A lU 


1 . UH- A 1 U 


-1 


f\ 1 n V in 


-2 




2.00 X 10'^^ 


3.23 X 10- 


-3 


1.04 X 10- 


-3 


ISO (NED) 


4.44 X 10'^^ 


7.90 X 10- 


-4 


2.60 X 10- 


-4 


ISO (NED) 


3.33 X 10'* 


5.00 X 10- 


-5 






(Farrah et aIJ2002<0 


3.28 X 10'* 


7.97 X 10- 


-5 


4.34 X 10- 


-6 


SDSS-DR5 


3.93 X 10'* 


9.09 X 10 


-5 


1.26 X 10- 


-6 


SDSS-DR5 


4.28 X 10'* 


1.14 X 10 


-4 


1.05 X 10- 


-5 


APM 


4.81 X 10'* 


1.01 X 10- 


-4 


1.02 X 10- 


-6 


SDSS-DR5 


6.28 X 10'* 


8.52 X 10- 


-5 


7.85 X 10- 


-7 


SDSS-DR5 


8.47 X 10'* 


8.47 X 10- 


-5 


2.26 X 10- 


-6 


SDSS-DR5 



Table B.ll. Photometric data for IRAS F15307+3252. 



y (Hz) 




f V (Jy) 




Error (Jy) 




Ref. 




1.40 X 


10" 


5.71 X 


10 


-3 


1.09 X 10- 


-4 


VLA (NED) 




8.42 X 


10' 


9.20 X 


10 


-4 


4.00 X 10- 


-5 


VLA (NED) 




1.02 X 


10" 


4.50 X 


10- 


-2 






OVMA (NED) 




2.39 X 


10" 


5.10X 


10- 


-3 






OVMA (NED) 




3.53 X 


10" 


1.15 X 


10 


-2 






SCUBA (ROTah et d:2()02a|) 


6.66 X 


10" 


1.06 X 


10 


-1 






SCUBA (Farrah et aIJ2002d) 


1.66 X 


10'2 


4.14 X 


10- 


-1 


1.76 X 10- 


-1 


ISO (NED) 




3.00 X 


10'2 


5.10X 


10 


-1 


6.20 X 10- 


-2 


IRAS (NED) 




3.15 X 


10'^ 


3.68 X 


10 


-1 


1.16 X 10- 


-1 


ISO (NED) 




5.00 X 


10'2 


2.80 X 


10- 


-1 


2.70 X 10- 


-2 


IRAS (NED) 




1.20 X 


10'^^ 


8.00 X 


10- 


-2 


2.40 X 10- 


-2 


IRAS (NED) 




1.26 X 


10'^' 


5.70 X 


10 


-2 


3.00 X 10- 


-3 


Spitzer -IRAC 




2.50 X 


10'3 


4.50 X 


10 


-2 






IRAS (NED) 




1.39 X 


10'* 


3.20 X 


10- 


-4 


3.00 X 10- 


-5 


MMT (NED) 




1.80 X 


10'* 


1.68 X 


10 


-4 


1.40X 10- 


-5 


MMT (NED) 




2.43 X 


10'* 


1.48 X 


10 


-4 


1.20 X 10- 


-5 


MMT (NED) 




3.33 X 


10'* 


6.70 X 


10- 


-5 


1.85 X 10- 


-5 


SSS 




3.28 X 


10'* 


1.28 X 


10- 


-4 


3.91 X 10- 


-6 


SDSS-DR5 




3.93 X 


10'* 


7.17 X 


10 


-5 


1.06 X 10- 


-6 


SDSS-DR5 




4.28 X 


10'* 


6.32 X 


10 


-5 


1.75 X 10- 


-5 


SSS 




4.81 X 


10'* 


5.90 X 


10- 


-5 


8.70 X 10- 


-7 


SDSS-DR5 




6.28 X 


10'* 


4.40 X 


10 


-5 


6.89 X 10- 


-7 


SDSS-DR5 




6.81 X 


10'* 


3.40 X 


10 


-5 


4.00 X 10- 


-6 


XMM-Newton 


-OM 


8.21 X 


10'* 


2.90 X 


10- 


-5 


2.00 X 10- 


-6 


XMM'Newton 


-OM 


8.47 X 


10'* 


3.34 X 


10- 


-5 


1.69 X 10- 


-6 


SDSS-DR5 




9.67 X 


10'* 


3.00 X 


10 


-5 


2.00 X 10- 


-6 


XMM-A'evvton 


-OM 


1.48 X 


10'5 


1.90 X 


10- 


-5 


6.00 X 10- 


-6 


XMM-Newton 


-OM 



A. Ruiz et al.: Spectral Energy Distribution of Hyper-Luminous Infrared Galaxies, Online Material p 4 



Table B.12. Photometric data for IRAS 16347+7037. Table B.13. Photometric data for IRAS 18216+6418. 
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